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This paper studies short-range order (SRO) in the semiconductor alloy (GaN)i_a;(ZnO)j;. Monte 
Carlo simulations performed on a density functional theory (DFT)-based cluster expansion model 
show that the heterovalent alloys exhibit strong SRO because of the energetic preference for the 
valence-matched nearest-neighbor Ga-N and Zn-O pairs. To represent the SRO-related structural 
correlations, we introduce the concept of Special Quasi-ordered Structure (SQoS). Subsequent DFT 
calculations reveal dramatic influence of SRO on the atomic, electronic and vibrational properties of 
the (GaN)i_ 2 :(ZnO) 2 : alloy. Due to the enhanced statistical presence of the energetically unfavored 
Zn-N bonds with the strong Zn3d-N2p repulsion, the disordered alloys exhibit much larger lattice 
bowing and band-gap reduction than those of the short-range ordered alloys. Inclusion of lattice 
vibrations stabilizes the disordered alloy. 
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I. INTRODUCTION 


Relative to pure end-member materials, the non-isovalent pseudobinary semiconductor alloy (GaN)i_ 2 ,(ZnO)a; shows 
improved efficiency as a photocatalyst in splitting water into hydrogen and oxygen under visible light illuminatiorf^. 
High efficiency is partly attributed to the band-gap reduction which can be tuned by varying the ZnO content x 
of the alloy. First-principles calculations show that including short-range order (SRO) affects the band gaps of the 
heterovalent semiconductor alloyd^is] Experiments on different (GaN)i_a;(ZnO)a; samples also observe large variation 
in the band gaps, which could be attributed to the different degrees of SRO introduced in growing the samples. For 
example, the absorption edge shifts monotonically to longer wavelength with increasing x for samples synthesized by 
nitridation of nanocrystalline ZnGa 204 and ZnO precursor^, while a minimum gap at a; ^ 0.5 is found for samples 
synthesized by mixing of GaN and ZnO powders at high pressure and high temperatur^. Despite the experimental 
indication of the presence of SRO, a thorough theoretical understanding is still lacking. An even m ore challenging 
question is how the vibrational properties depend on SRO, and how they influence the degree of SRd^. The effect 
of lattice vibrations is important for isovalent semiconductor alloy thermodynamic calculation^. To our knowledge, 
to date, there are no published phonon data for (GaN)i_a;(ZnO)a;. Therefore we perform detailed first-principles 
investigations to illuminate the role SRO plays on phonons in the (GaN)i_a;(ZnO)a; alloy. 

Previous theoretical studie s assu me the (GaN)i_a;(ZnO)a, alloy to be completely randorcPtll] The Special Quasiran¬ 
dom Structure (SQS) methoc^^ is often used to construct supercells mimicking random alloy^l. However, even for 
isovalent ternary semiconductor alloys, neglecting SRO introduces non-negligible systematic error^l^lll^. For example, 
the band gaps of Alo.sGao.sAs, Gap.sing, 5 P and Alo. 5 Ino. 5 As alloys can be reduced by as much as O.leV through 
clusterin j^^b7 | ^pEe electronic properties of wurtzitic Gai_a;Ina;N and Ali_a;Ina;N alloys are also found to be very 
sensitive to SRO in the cation distributioiP^. For quaternary alloys the sensitivity to SRO is even greateJ^SEni. The 
situation is compounded for the (GaN)i_a;(ZnO)a; alloy whose heterovalent nature favors local charge neutrality and 
therefore valence-matched nearest-neighbor Ga-N and Zn-0 pairs. In our previous study on the (GaN)i_a,(ZnO)a; 
alloji^, referred to hereafter as I, first-principles calculations combined with the cluster expansion methocP^HlIl ^nd 
Monte Carlo simulations predicted a large degree of SRO. In constructing the cluster expansion model, the total 
energy of a specific configuration is calculated in its relaxed structure. Local relaxations of surprisingly large mag¬ 
nitude are found in our subsequent studji^. The aim of the present study is to construct DFT-affordable supercells 
whose structural correlations accurately reflect the SRO found by the above approach. The method we use (“special 
quasi-ordered structure, SQoS”) was used in 1998 by Saitta et aP^but rarely since theiP^. This method allows us to 
study with a single DFT calculation, for each x, the statistical average atomic, electronic and vibrational properties of 
the (GaN)i_a;(ZnO)a; alloy. The dependence of structural properties such as bond-length distribution and bond-angle 
variation upon SRO will be discussed in a separate paper^Zl. 


II. COMPUTATIONAL METHOD 


The (GaN)i_a;(ZnO)a; alloy was modeled in wurtzite structure with interpenetrating cation and anion hep sub¬ 
lattices. Ga/Zn can only occupy the cation sublattice, while N/0 can only occupy the anion sublattice. A detailed 
description of the cluster expansion model used in this study can be found in I. Monte Carlo simulations are performed 
using the AT AT package^Sn^ with a 12 x 12 x 8 supercell containing 4608 atoms. For each (a;,T), an ensemble of N 
configurations (labeled by s = 1,2, ...,N) is equilibrated for I x 10^ MC passes followed by a subsequent 1 x 10"* MG 
passes sampling. The site occupation is denoted by Ising spin ai with cr = I denoting Ga/N and a = —1 denoting 
Zn/0 respectively. For the structural correlations, the notations are adopted from Ref. [T^ The total energy of 
sample s is expanded in terms of clusters (called “figures” and labeled as {k, m}). The label k = 1, 2,... is the number 
of sites of the cluster. The label m = 1, 2,... enumerates the distinct cluster geometries, as shown in the inset of Fig. 
1. The structural correlation function describes the occupation of the cluster. The label (Z,s) indicates 

that the cluster is located at location I in sample s. For instance, the value of n 2 p(Z,s) is the spin product aiaj 
for the particular nearest-neighbor pair of sites i and j positioned at location I in sample s. nfc^m(s) refers to the 
average of Hj. ,„(/, s) over all locations I in one sample s, and (H/j „,) refers to the average of H^, „,(s) over the samples 
equilibrated at a certain {x,T). 

The motivation of the SQS approadP^EH is to approximate the actual alloy with one representative special structure 
S whose struct ural correlation functions nfc_m(S') best match the corresponding ensemble-averaged (Hfc^m) of the real¬ 
istic allojffimi The origi nal S QS approach reproduces the average structural correlation functions of the random (i?) 
alloy nfc,m(<S') ~ We extend the SQS approach to the correlation functions of short-range ordered alloys 

{Ak,m)sRO- We first obtain by performing Monte Carlo simulations on a DFT-based cluster expansion 

model. Then we generate numerous site occupancies for a certain composition x and look for the representative config- 
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uration S for which the set of nfe_m(-S') is closest to (nk,m)sRO by minimizing J2k,m 9k,mDk,m [nfc,m(s) - 
where Dk^m is the degeneracy (number of equivalent figures) and gk,m is the assigned weighting factor. Enumeration 
of all possible configurations is not possible since the number grows exponentially with the number of atoms in the 
supercell. However, increasing the size of the supercell allows better flexibility of matching structural correlation 
functions. The conflict is eased by the short-range nature of the structural correlations of the (GaN)i_a;(ZnO)a; alloy. 
The most relevant physical property of the (GaN)i_a;(ZnO)a; alloy is the formation energy, which is dominated by 
the short-range pair structural correlation^^. We thus assign large weighting factors to the nearest-neighbor {2,1-2} 
(meaning {2,1} and {2,2}) and next nearest-neighbor {2,3-6} figures. At each {x,T) we generate 1 x 10® 72-atom 
supercell (3x3x2) candidate structures among which the best-matching structure S is chosen. We emphasize that 
the obtained structures are not the optimal SQS. However, the contributions to the energetics from longer-range 
figures E{s) - (E) = Y.k,m^k,m [nfc,m(s) - (nfe,m)sijo] e/c,m (Eq. (3.1) in Ref. HSl) are reasonably small. To avoid 
confusion, we name the corresponding special structure for the short-range ordered alloy (equilibrated at the ex¬ 
perimental synthesis temperature T = 1, 123KP) as SQoS (Special Quasi-ordered Structure), and for the disordered 
alloy (equilibrated at an unrealistic high temperature T = 20,000K) as SQdS (Special Quasi-disordered St ructu re), 
in resemblance to the widely used SQS (Special Quasi-random Structure) formalism introduced by Zungeil^^ESl, A 
completely random “SQS” at x = 0.5 is also studied for reference. The constructed special structures are provided in 
supplementary materiald^. 

The constructed special structures are fully relaxed with respect to atomic coordinates, volume and shape. Elec¬ 
tronic structure calculations are performed using the Quantum ESPRESSO packag^^with the PBE sol functionaP^. 
The pseudopotentials are constructed by means of the projector augmented wave (PAW) methocP^*^ with 60 Ry and 
240 Ry cutoff energy for plane-wave basis set and charge density respectively. Ga-3d and Zn-3fi states are treated 
explicitly as valence states. The A:-point mesh is chosen to be equivalent to a 6 x 6 x 4 mesh for the 4-atom wurtzite unit 
cell. To speed the structural relaxations, the input lattice parameters are estimated using Vegard’s lavP^. Nowadays 
fairly large supercells (e.g., over 50 atoms) can be handled at the DFT level. Eor the non-isovalent semiconductor 
alloys where large structural relaxations are expected, one can benefit greatly in terms of the computational efficiency 
from a pre-relaxation prior to the expensive DFT total energy and force calculations. We will address the issue of 
pre-relaxation in a subsequent stud}P^. Phonons are calculated using the small displacement method as implemented 
in the PHON cod^^. For each 72-atom primitive cell, a 2 x 2 x 2 supercell is used while a small displacement of 0.02 
A is employed. The force constants are calculated with the SIESTA packag^^. Pseudopotentials for all the atomic 
species are available from the SIESTA homepage http://www.icmab.es/siesta/^, except for Ga a smaller d-orbital 
cutoff radius is usecP^. 

GaN and ZnO have a type-H band alignmenflii]. The valence band is composed mainly of N-2p states. DFT with 
LDA or GGA tends to over-delocalize the semicore Zn-d states and consequently over-hybridize the semicore Zn-d 
states with the N-p states, resulting in an enhancement of the p-d repulsion. The band gap is therefore severely 
underestimated due to the artificially large p-d repulsion. In this study we add U corrections to the semicore Ga-d 
and Zn-d state^^. The on-site Coulomb interaction parameter U ~ 3.1eV is determined by a first-principles method 
adopted in Ref. [42l U is approximated as the screened atomic on-site Coulomb interaction [/“‘/coo, where C/“‘ is the 
Coulomb energy cost of placing two electrons at the same site (C/“* = Eat{d"'~^^) -P Eat{d"'~^) — 2Eat{d^)) and too is 
the optical (high-frequency) dielectric constant. In this study we take d® occupancy as the reference point for d” and 
evaluate from DFT atomic energies. The optical dielectric constant Coo is calculated from linear-response theorj^^. 
A similar approach of screening the exact-exchange by the dielectric constant is shown to improve significantly the 
performance of the traditional hybrid functionals^. 


III. RESULTS AND DISCUSSIONS 

A. Structural correlation 

As is predicted in I, the (GaN)i_a;(ZnO)a; alloy is thermodynamically stable over the full range of compositions 
for T > 870K. The x = 0.5 alloy orders at low temperature and undergoes a first-order order-disorder transition at 
T « 870K. The ground state is an ordered 50%-50% superlattice labeled as (GaN)i(ZnO)i, where GaN and ZnO 
double layers stack alternately along the hexagonal c-axis (PG^mc). The formation energy for the (GaN)i (ZnO) i 
superlattice is predicted to be small and negative, indicating weak stability against phase separation. An analogous 
superlattice structure is also predicted for the (SiC)m(AlN)„ alloji^. 

Upon alloying, the main effect of SRO is to enhance the statistical presence of the valence-matched nearest-neighbor 
Ga-N and Zn-0 pairs. The ensemble-averaged pair correlation functions sd x = 0.5 (Fig. 1) reveal a large degree 

of SRO. The nearest-neighbor (n 2 ,i_ 2 ) deviate significantly from the null value of the random alloy, while the next 
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nearest-neighbor ( 112 ,s-e) are relatively small, comparable with those found in ternary nitride isovalent semiconductor 
alloyJi^. Longer-range ( 112 , 7 - 14 ) are not important. The long tail of the (nfc,^) — T curve also indicates that SRO 
persists to high temperature, and therefore complete randomness may not be achievable under common experimental 
growth conditions. The positive signs of (112, 1 - 2 ) indicate nearest-neighbor preference for the valence-matched Ga-N 
and Zn-0 pairs, while the positive signs of ( 112 ,a-e) indicate next nearest-neighbor preference for Ga-Ga and Zn-Zn 
as well as N-N and 0-0 pairs. The composition dependence of (11^,^) at T = 1,123K is shown in Fig. 2. The 
deviation of (nfe,^) from that of the random alloy increases upon mixing, and yields the largest deviation at x = 0.5, 
where neglect of SRO is worst. To compare the degree of SRO included in SQoS, SQdS and SQS, we summarize in 
Table 1 the corresponding structural correlation functions at x = 0.5. The 72-atom SQoS, SQdS and SQS accurately 
reproduce the ensemble-averaged structural correlation functions obtained with a 12 x 12 x 8 supercell. These special 
structures are expected to yield an accurate description of the atomic, electronic and vibrational properties of the 
(GaN)i_a;(ZnO)a; alloy. 



Temperature (K) 

FIG. 1. Ensemble-averaged pair correlation functions (112,m) at x = 0.5. Definitions of pair figures {2,m} can be found in I. 
{2,1-2} and {2, 3-6} stand for nearest-neighbor cation-anion pair figures and next nearest-neighbor cation-cation/anion-anion 
pair figures respectively. Longer-range pair figures {2, 7-14} are shown by dotted lines. 



ZnO content x 

FIG. 2. Ensemble-averaged pair correlation functions (112,™) at T = 1,123K. The structural correlations for the random alloy 
(nfc,™)_R = (2x — 1)*^ is shown by the solid grey line for comparison. 
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TABLE I. 72-atom SQS at a; = 0.5. In spite of the small size of the supercell, the statistical accuracy is good. 



Ba.rn-SQOS 

{n.2,m)l,123K 

Ha.^-SQdS 

(n2,m)20,000X 

n2,,„-SQS 

{2,1} 

0.444 

0.442 

0.074 

0.070 

0 

{2,2} 

0.333 

0.333 

0.000 

0.058 

0 

{2,3} 

0.037 

0.041 

-0.037 

-0.011 

0 

{2,4} 

0.037 

0.036 

0.000 

-0.012 

0 

{2,5} 

0.074 

0.089 

-0.037 

-0.012 

0 

{2,6} 

0.074 

0.063 

0.000 

-0.016 

0 


TABLE 11. Calculated e^o and the corresponding U parameters for GaN and ZnO. Experimental values are shown in 
parenthesis. The PBE versiorP^of the GGA functional is used instead of PBEsol in obtaining t/“*, due to its better treatment 
of free atoms. 



[/“* (eV) 

€oo 

U (eV) 

GaN 

18.1 

5.9 (5.35) 

3.1 

ZnO 

16.1 

5.2 (3.71) 

3.1 


B. Atomic, Electronic and Vibrational Properties 

The calculated 17“*, e^o and U parameters are listed in Table 2. Compared to the experimental value^^, the 
calculated optical dielectric constant is overestimated due to the band-gap underestimation of DFT. However, since 
the atomic and electronic structures of GaN and ZnO are not very sensitive to the U parameters, the error in the 
calculated e^o (and also the choice of the reference point for d”) does not affect the main conclusions drawn in this 
study. The calculated lattice constants and band gaps are listed in Table 3. DFT-PBEsol calculations accurately 
reproduce the lattice constants of GaN and ZnO. The band gap of ZnO is more sensitive to the U correction, due to 
the strong interaction between the high-lying Zn-Sd states and the 0-2p states. We then perform DFT-l-17 calculations 
on the SQoS, SQdS and SQS in order to obtain accurate electronic structure properties. For comparison, total energy 
and force calculations on configurations randomly selected from the T = 1 ,12377 ensembles are also performed within 
the DFT+C/ methodology. As shown in Fig. 3, the constructed SQoS accurately represents the ensemble-averaged 
energetics of the short-range ordered (GaN)i_a;(ZnO)a; alloy. The formation energy of SQoS is significantly lower than 
that of SQdS. The effect of SRO on the energetics grows upon mixing. Even at T = 20, OOOK, the formation energy 
of SQdS is still considerably lower than that of SQS due to the non-negligible residual SRO. 

SRO also plays an important role in determining the structural properties. Fig. 4 compares the lattice constant 
bowing obtained theoretically and experimentally. Once again, the lattice constants of SQoS accurately reproduce 
the corresponding ensemble-averaged values. With reduced SRO, the disordered alloy shows an expa nsio n as well as 
a larger bowing compared to the short-range ordered alloy. The experimentally synthesized sampled^ also exhibit 
moderate bowing, larger than the short-range ordered alloy but smaller than the disordered alloy, indicating the 
presence of SRO. 

Fig. 5 shows the (nearest-neighbor) bond-length distribution of the short-range ordered (T=1123K) alloy at x = 0.5. 
In the (GaN)i_a;(ZnO)a; alloy, the Ga-N bonds shrink while the Zn-0 bonds expand. This unusual bond-length 
distribution is determined by the non-isovalent nature of the alloy. A follow-up studj^^ will discuss the prediction 
and explanation of the bond-length distribution based on the concept of bond valenc^^. For the (GaN)i_ 2 ,(ZnO)a; 
alloy, the Zn-N bond-length distribution has crucial importance since it is related to the band-gap reduction through 



TABLE III. 

Calculated lattice constants 

a and 

c and band gaps Eg 

for GaN and ZnO. 




GaN 



ZnO 



a (A) 

c(A) 

E, 

a (A) 

c (A) 

Eg (eV) 

PBEsol 

3.182 

5.187 

1.88 

3.225 

5.207 

0.71 

PBEsol-tt/ 

3.184 

5.189 

1.89 

3.232 

5.213 

0.92 

Expt. 

3.189 

5.185 

3.3 

3.250 

5.204 

3.4 
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FIG. 3. DFT-calculated formation energies of SQoS, SQdS and SQS. {SQoS} is a set of configurations randomly selected from 
the T — 1,123K ensemble. 


the Zn3fi-N2p repulsion. In Fig. 6 we show its dependence on the ZnO content. We find that the Zn-N bond-length 
distribution of the short-range ordered alloy shifts to shorter bonds as the ZnO content increases. Shorter Zn-N bond 
lengths result in stronger Zn3d-N2p repulsion and therefore significantly push up the top of the valence band. 

The bond-angle variation is also unusual, namely N-Ga-N and Ga-N-Ga angles expand while O-Zn-0 and Zn-O-Zn 
angles shrink relative to the ideal tetrahedral angle 109.5°. Fig. 7 shows the variation of bond angles. For example, 
the Ga centered bond angle shrinks with increased presence of ligand O atoms. This tendency can also be explained 
using the concept of bond valence. For Fig. 5-7, see Ref. [27] for a statistically reliable prediction based on the bond 
valence method. 



o 

c 

V) 

C 

O 

o 

(D 

O 


FIG. 4. DFT-calculated lattice constants of SQoS and SQdS. 


The atomic and electronic structures of the short-range ordered alloys deviate significantly from those of the 
disordered alloys. Theoretical atomistic modeling requires explicit inclusion of SRO. Fig. 8 shows the calculated band 
gaps of SQoS and SQdS. Due to the enhanced statistical presence of the Zn-N bonds, the band gap of the disordered 
alloy is further reduced relative to that of the short-range ordered alloy. The band-gap reduction is asymmetric. For 
the disordered alloys the band gap bowing is somewhat parabolic, while for the short-range ordered alloys the band 
gap reduces almost linearly with increasing ZnO in the GaN host. The linear band-gap reduction is maintained even 
for the unrelaxed short-range ordered alloys, indicating the dominating role of configurational SRO. In Fig. 8 we 
also show the linear redshift of the absorption onset with increased ZnO content observed in samples synthesized 
by nitridation of nano crystalline ZnGa 204 and ZnO precursor^ The linearity is a clear indication of the presence 
of SRO. We also notice that the high-temperature and high-pressure synthesized samples exhibit the minimum gap 
at a: = 0. P, which is somewhat consistent with the parabolic band gap bowing of the random alloy. The parabolic 
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* ' - L*^ • A X — — I N, - - _ 

1.8 1.9 2.0 2.1 2.2 

Bond length (A) 


FIG. 5. DFT-calculated bond-length distribution at T = 1123K and x = 0.5. Thirty 72-atoni structures are selected from 
the corresponding thermodynamic ensemble. The bin interval is set to O.OlA. The vertical lines mark the bond lengths of the 
corresponding compounds. 



- . _I_,_ I I , -- 

1.85 1.90 1.95 2.00 2.05 2.10 

Bond length (A) 


FIG. 6. DFT-calculated Zn-N bond-length distribution a.t T = 1123K. 


bowing is attributed to the promoted kinetics of mixing at high-temperature and high-pressure. The contrast in the 
band gap bowing is a clear indication of the importance of SRO. Since the SRO introduced in the sample is related 
to the synthesis techniques and the growth conditions, one might therefore consider the opportunity of engineering 
the band gap Eg{x,T,Ii) via SRO. 

Fig. 9 compares the projected density of states (PDOS) of SQoS and SQdS at a; = 0.5. The main contribution 
to the bottom of the valence band comes from 0-2s states, which are taken as the reference level because they are 
less sensitive to the local chemical environment. The top of the valence band is mainly composed of N-2p states. 
For the disordered alloy the increased statistical presence of the energetically unfavored Zn-N pairs pushes the band 
edge upward, resulting in further reduction of the band gap. The N-2p states depend strongly on the local chemical 
environment. Fig. 10 shows the PDOS of N-2p states with the N atoms surrounded by different numbers of Zn atoms. 
The N-2p states shift upward (dashed lines in Fig. 10) with increased presence of Zn neighbors. 

The effect of lattice vibrations is calculated using the harmonic approximation. The phonon DOS for the SQoS 
alloys along with those of compound GaN and ZnO are shown in Fig. 11. Three mechanisms have been suggested to 
explain the origin of vibrational entropy differences in alloy^^ the bond proportion effect, the volume effect and the 
size mismatch effect. Upon disordering, the bond proportion effect is typically associated with a broadening of the 
phonon DOS due to the statistical presence of bonds with different stiffness. On the other hand, the volume effect is 
usually characterized by an overall shift of the phonon DOS due to the change in the frequency of all phonon modes. 
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Bond angle 0 (°) 

FIG. 7. DFT-calculated bond-angle variation at (T = 1123K,2: = 0.5). 



ZnO content x 


FIG. 8. DFT-calculated band gaps of SQoS and SQdS. Experimental measurements (Lee et al in Ref. |4]and Chen et al in Ref. 
[S]) are also shown for comparison. 


For the (GaN)i_ 2 ,(ZnO)a; alloy, as the alloy expands with increasing x, the phonon DOS shifts to lower frequencies as 
the chemical bonds are in general softened. The volume effect is magnified by the fact that the “ionic” Zn-0 bond is 
softer than the “covalent” Ga-N bond. The low-frequency phonon DOS is well represented by the composition weighted 
average (1 — x)gGaJ^iuj)+xgznoi^)- A significant part of the phonon DOS difference (and therefore the vibrational 
entropy difference) comes from the high-frequency phonons. The effect of SRO is shown for the x = 0.5 case. The 
high-frequency phonon DOS of the SQdS exhibits a much broader spectrum than that of the SQoS. Gonsequently the 
phonon mixing entropy of the SQdS is three times larger than that of the SQoS, as is shown in Fig. 12. While the 
cc-dependence of the configurational mixing entropy is symmetricP^, the cc-dependence of the phonon mixing entropy 
is highly asymmetric, indicating that the inclusion of the vibrational free energy into the alloy thermodynamics could 
alter the shape of the phase diagram. 


IV. CONCLUSIONS 

The importance of SRO in atomistic modeling schemes such as the SQS approach is often overlooked. For bi¬ 
nary metal alloys or isovalent semiconductor alloys, SRO is usually less important. However, for the non-isovalent 
semiconductor alloys, the valence-matching driving force induces signihcant SRO. The SQS approach provides a way 
of approximating the actual alloy with a DFT-affordable supercell. In order to properly compute the non-isovalent 
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Energy (eV) 


FIG. 9. Projected density of states (PDOS) of the valence band. The cation-s states in the conduction band are not shown. 
The deep-lying anion-s states are shown by the shaded area. 



Energy (eV) 


FIG. 10. PDOS of N-2p with N atoms being surrounded by 0, 1, 2, 3 or 4 Zn neighbors. PDOS is not normalized so that the 
area under each curve represents the corresponding statistical presence. 


alloy, one needs prior knowledges of SRO. In this study the correlated site occupations are provided by Monte Carlo 
simulations on a DFT-based cluster expansion model. Exhaustive enumeration of all site occupations is avoided due 
to the SRO in the (GaN)i_ 2 ,(ZnO)a; alloy. We seek to match only the cation-anion nearest-neighbor and the cation- 
cation/anion-anion next nearest-neighbor correlations. The longer-range correlations are optimized to a lesser extent. 
Since the short-range nature is generic in the non-isovalency, the construction of SQoS proposed in the present study 
should also be applicable to other non-isovalent semiconductor alloys. If longer-range correlations come into play, one 
might apply for example the evolutionary algorithm in order to efficiently search for the optimal SQoS. 

The present study reveals the presence of strong SRO in the (GaN)i_a;(ZnO)a; alloy. We construct reliable SQoS 
and SQdS whose structural correlations reproduce those of the short-range ordered alloys and the disordered alloys 
respectively. Atomic, electronic and vibrational properties of the short-range ordered alloys deviate significantly from 
those of the disordered alloys. The short-range ordered alloys experience smaller lattice bowing than the disordered 
alloys. We offer a tentative explanation in terms of SRO for the discrepancy of the band gaps found in samples 
synthesized by different methods. SRO inhibits the nearest-neighbor Zn-N pairs, which affects the strength of the 
Zn3d-N2p repulsion and consequently the top of the valence band. The dependence of the N-2p states on local chemical 
environment demonstrates the vital role of SRO in accurately describing the (GaN)i_ 2 ;(ZnO)a; alloy. Phonon DOS is 
sensitive to the presence of SRO. Disordered alloys have much larger vibrational entropy of mixing than short-range 
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Frequency (cm’’) 


FIG. 11. Phonon DOS for the SQoS (solid black) {x = 0.1,0.3, 0.5, 0.7 and 0.9) and SQdS (dash-dot blue) (x = 0.5) alloys. 
The dash red lines represent the corresponding average of phonon DOS: (1 — x)gGaN(‘jj)+xgzno(‘^)- 



FIG. 12. x-dependence of the phonon mixing entropy ASvib = Svibix) — [(1 — x)S'® 5 ’^ -I- xS'^f^]. 


ordered alloys. 
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